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Chiral symmetry restoration at nonzero temperature is studied in the framework of the 0(4) linear 
sigma model and the U(Nf) r x U(Nf)i linear sigma model with Nf = 2, 3 and 4 quark flavors. We 
investigate the temperature dependence of the masses of the scalar and pseudoscalar mesons, and 
the non-strange, strange, and charm condensates within the Hartree approximation as derived from 
the Cornwall- Jackiw-Tomboulis formalism. We find that the masses of the non-strange and strange 
iy-^ ' mesons at nonzero temperature depend sensitively on the particular symmetry of the model and the 

qq I number of light quark flavors Nf. On the other hand, due to the large charm quark mass, neither do 

charmed mesons significantly affect the properties of the other mesons, nor do their masses change 
appreciably in the temperature range around the chiral symmetry restoration temperature. In the 
chiral limit, the transition temperatures for chiral symmetry restoration are surprisingly close to 
£C) , those found in lattice QCD. 

o 



q ■ I. INTRODUCTION 

For Nf massless quark flavors, the QCD Lagrangian has a chiral U(Nf) r x U(Nf)e = SU(Nf) r x 
SU(Nf)e x U(l)v x U(1)a symmetry. Here, V = r + £, while A = r — £. The U(l) v symmetry 
corresponds to baryon number conservation. It is always respected and thus plays no role in the 
r> I symmetry breaking patterns considered in the following. In the vacuum, a non- vanishing expectation 

value of the quark condensate, (qi q r ) ^ 0, spontaneously breaks the above symmetry to SU(Nf)v- 
This gives rise to Nf Goldstone bosons which dominate the low-energy dynamics of the theory. As 
shown by 't Hooft 0,0, instantons break the U(1)a symmetry explicitly to Z{Nf)A (For the 
low-energy dynamics of QCD, however, this discrete symmetry is irrelevant.) Consequently, one 
of the N"j Goldstone bosons becomes massive, leaving N"j — 1 Goldstone bosons. The SU(Nf) r x 
SU(Nf)e x U(1)a symmetry of the QCD Lagrangian is also explicitly broken by nonzero quark 
masses. The N* — 1 low-energy degrees of freedom then become pseudo-Goldstone bosons. For 
M < Nf degenerate quark flavors, an SU(M)v symmetry is preserved. 

As indicated by lattice QCD calculations y|, chiral symmetry is restored at temperatures around 
~ 150 MeV (at zero net- baryon number density). It is difficult to determine the order of the chiral 
phase transition on the lattice. At this time, lattice calculations have not unambiguously answered 
this question. For physical values of the quark masses, calculations with staggered fermions [5J favor 
a smooth crossover transition, while calculations with Wilson fermions Q predict the transition to 
be of first order. 

For vanishing quark masses, i.e., in the chiral limit, however, one can use universality arguments 
to determine the order of the phase transition. According to universality, the order of the chiral 
transition in QCD is identical to that in a theory with the same chiral symmetries as QCD, for 
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instance, the U(Nf) r x U(Nf)g linear sigma model. This argument was employed by Pisarski and 
Wilczek Q who found that for Nf = 2 flavors of massless quarks, the transition can be of second or- 
der, if the U symmetry is explicitly broken by instantons. It is driven first order by fluctuations, 
if the U(1)a symmetry is restored at T c . For Nf — 3 massless flavors, the transition is always first 
order. In this case, the term which breaks the U(1)a symmetry explicitly is a cubic invariant, and 
consequently drives the transition first order. In the absence of explicit U(1)a symmetry breaking, 
the transition is fluctuation-induced of first order. For Nf = 4, the same argument leads to a first 
order chiral transition in the absence of the U(1)a anomaly. The term which breaks the U(1)a 
symmetry is no longer a cubic invariant, but of quartic order in the fields. Since this term does not 
generate an infrared-stable fix point, the transition remains of first order. 

For nonzero quark masses, the chiral symmetry of QCD is explicitly broken. Nonzero quark 
masses act like a magnetic field in spin systems, such that a second order phase transition becomes 
a crossover transition. When the quark masses increase, a first order phase transition may for 
a while remain of first order, but it will ultimately become a crossover transition, too. In order 
to decide whether this happens for a particular choice of quark masses, universality arguments 
cannot be applied and one has to resort to numerical calculations. As an alternative to lattice QCD 
calculations, one can also use linear sigma models to investigate this question 0. 

Studying the linear sigma model at nonzero temperature, however, requires many-body resum- 
mation schemes, because infrared divergences cause naive perturbation theory to break down |8(]. In 
these resummation schemes, one necessarily has to make approximations by selecting certain subsets 
of diagrams. The most commonly employed scheme is the Hartree approximation. This approxima- 
tion fails to reproduce the correct order of the chiral transition for the 0(4) linear sigma model: the 
transition is of second order, while the Hartree approximation yields a first order transition. For 
U(Nf) r x U(Nf)e linear sigma models with Nf = 3, however, the Hartree approximation correctly 
produces a first order chiral transition [jj. For Nf = 2, in the absence of the U(1)a anomaly the 
transition is of first order and, as will be shown in this paper, the Hartree approximation agrees 
with this result. With the U(1)a anomaly, the transition is of second order. We shall find that the 
Hartree approximation (slightly) violates this prediction by producing a (weak) first order transition. 

The utility of linear sigma models, however, transcends further than just predicting the order of 
the chiral phase transition. The degrees of freedom of the U(Nf) r x U(Nf)g linear sigma model are 
the scalar and pseudoscalar mesons. In the vacuum, the latter are the (pseudo-) Goldstone bosons 
of chiral symmetry breaking. Thus, the linear sigma model can also be viewed as an effective low- 
energy theory for QCD. At high temperatures, chiral symmetry is restored. Consequently, chiral 
partners among the scalar and pseudoscalar mesons must become degenerate in mass. Since the 
U(Nf) r x U(Nf)i linear sigma model treats both scalar and pseudoscalar mesons on the same footing, 
it is particularly suited to describe the change of meson properties across the chiral transition. Note 
that, in the case of a first order chiral transition which coincides with the deconfinement transition, 
the correct degrees of freedom in the high-temperature phase are quark and gluons instead of mesons. 
However, in the case of a crossover transition mesonic degrees of freedom are well-defined even above 
the chiral transition temperature. 

In the present work, we follow this line of arguments and investigate the change of meson masses 
and quark condensates with temperature in the framework of U (Nf) r x U (Nf)g linear sigma models. 
We focus in particular on the question how the number of quark flavors Nf changes the temperature 
dependence of the meson masses and the quark condensates. For a given Nf, we furthermore 
investigate the different patterns of symmetry breaking arising from the presence or absence of the 
U(1)a anomaly, and from taking zero or nonzero values for the quark masses. In this sense, the 
present study is an extension of previous work [fill Id). Table [I] presents an overview of the different 
models and patterns of symmetry breaking studied in this paper. 

The U(2) r x U(2)i linear sigma model has eight degrees of freedom: the scalar fields are the a 
meson and the three do mesons, the pseudoscalar fields are the r\ meson and the three pions. With 
spontaneous breaking of the U(2)a chiral symmetry, but without the U(1)a anomaly, the pions and 
the 1) meson are (pseudo-) Goldstone bosons, while the a and ao meson are heavy states. With the 
U(1)a anomaly, the rj meson also becomes heavy. 
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TABLE I: The table shows the symmetry breaking patterns studied in this paper. The y/ symbol indicates 
that the specific symmetry breaking pattern is physically interesting and investigated here. The x symbol 
indicates that this symmetry breaking pattern does not exist at all and the — symbolizes that the respective 
symmetry breaking pattern is not studied. 



With explicit breaking of the U(1)a symmetry due to the anomaly (and neglecting the U(l)v 
symmetry of baryon number conservation), the remaining symmetry of the U(2) r x U(2)i linear 
sigma model is SU(2) r x SU{2)^. This group is isomorphic to 0(4). The appropriate effective 
theory incorporating this symmetry is the 0(4) linear sigma model |ll| . This model has only four 
degrees of freedom, the a meson and the three pions. In this sense, it represents the limit of the 
U(2) r x U(2)i model for maximum U(1)a symmetry breaking. The U(2) r x U(2)t model is more 
general, since it allows to consider the properties of mesons also without the U(1)a anomaly. 

In the limit of maximum U(1)a symmetry breaking, the rj and ao mesons become infinitely heavy 
and are thus removed from the spectrum of physical excitations. This approximation is justified at 
small temperatures, where the dynamics is determined by the lightest hadronic degrees of freedom, 
i.e., the pions and, to some extent, the a meson. At higher temperatures and, in particular, around 
the chiral transition, however, heavier mesons become more and more important. It is therefore 
interesting to compare the results of the 0(4) model with those of the U(2) r x U(2)t model which 
incorporates these heavier degrees of freedom. Note that the 0(4) linear sigma model at nonzero 
temperature has been discussed in Refs. [hIIT^ . 

In the physical hadron spectrum, the kaons are lighter than the ao and r/ mesons and thus are more 
copiously produced at nonzero temperature. Consequently, the kaons are expected to influence the 
dynamics of the system to a larger extent than the ao and r/ mesons. Therefore, the U(2) r x U(2)i 
model is unrealistic in the sense that it neglects these strange meson degrees of freedom. From a 
physical point of view, it is thus necessary to enlarge the symmetry group to C/(3) r x [7(3)^ and 
study the corresponding linear sigma model. Nevertheless, it is still interesting to compare the 
U(2) r x U(2)i model with the f/(3) r x U{i)i model in order to see how the strange degrees of 
freedom affect the results. The U(3) r x U(3)i model was previously studied at nonzero temperature 
in Ref. 

Finally, we extend our investigations by including the charm degree of freedom. In principle, the 
corresponding linear sigma model has an SU(4) r x SU{4)i x U(1)a symmetry, but in nature this 
symmetry is strongly explicitly broken by the large charm quark mass. Therefore, we only consider 
the physically relevant case of explicit chiral symmetry breaking with U(1)a anomaly. 
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For all cases considered (see Table HJ we study the meson masses and the condensates as functions 
of temperature. In the chiral limit, we determine the phase transition temperatures by computing 
the effective potentials. In this aspect, we extend the previous study of Ref. jj. For the 0(4) model, 
this has already been done in Ref. [12 ■ 

Our calculations are done in the Hartree approximation which we derive within the Cornwall- 
Jackiw-Tomboulis (CJT) formalism The Hartree approximation has the advantage that the 

meson self-energies become independent of momentum and energy and one only has to solve gap 
equations for the meson masses as a function of temperature. To go beyond the Hartree approxi- 
mation, for instance by including energy-momentum dependent contributions to the self-energies, is 
considerably more difficult [Til [TBI fl6lll7| , and will be deferred to a future publication. 

The remainder of this paper is organized as follows. In Sec. |Hj we review the formulas of the 
CJT formalism, which are relevant for the present study, and apply them to the 0(4) model and 
the U(Nf) r x U{Nf)i linear sigma models for Nr = 2, 3, and 4. In Sec. II I II it is shown how 
to determine the coupling constants of the different models from the vacuum properties of the 
mesons and condensates. In Sec. IIVI we discuss the temperature dependence of the masses and the 
condensates. We conclude this work in Sec. Ivl with a summary of our results. 

We use the imaginary-time formalism to compute quantities at nonzero temperature. Our notation 

is 

jT/(A)=T £ J ^f(2mnT,k) , J J(x) = jf* * dr j 'd 3 x /(r,x) . (1) 

Our units are h = c = kg = 1. The metric tensor is g^ v = diag(+, — , — , — ). Throughout this work, 
all latin subscripts are adjoint U(Nf) indices, a = 0, . . . ,Nj — 1, and a summation over repeated 
indices is understood. 



II. LINEAR SIGMA MODELS IN THE HARTREE APPROXIMATION 



A. The CJT formalism 

The CJT formalism 0] generalizes the concept of the effective action for one-point functions to 
that of an effective action for one- and two-point functions. It is particularly useful for theories with 
spontaneously broken symmetry. In this case, truncating the standard loop expansion for the effec- 
tive action [l|| at some given order, at nonzero temperature one may obtain unphysical, tachyonic 
propagation of quasiparticles with small momenta. The reason for this failure of the standard loop 
expansion is that only the expectation value of the one-point function is self-consistently determined 
in this approach. The CJT formalism goes beyond the standard loop expansion by self-consistently 
determining the expectation value for the two-point function in addition to that of the one-point 
function. Effectively, this amounts to solving a Dyson-Schwinger equation for the two-point func- 
tions and yields a self-consistent computation of the quasiparticle self-energy. For translationally 
invariant systems, the effective action becomes the effective potential. To give an example, consider 
the general Lagrangian 

C{ip) = ^d^d^-U{ V ) (2) 
for a scalar quantum field ip. The effective potential in the CJT formalism reads 

V[4>,G] = U{4>) + \j\nG- 1 {k) + ^j [G-\k;$)G(k)-i\+VSA , (3) 

where is a c-number field (it is the expectation value of the quantum field ip in the presence of 
an external source), £/(</>) is the classical potential energy density (the tree-level potential) in the 
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Lagrangian J3J), and G 1 is the inverse of the tree-level propagator, 

G _1 (fc; 4>) = -k 2 + U"($) . 



(4) 



Here, U"((j)) is the second derivative of U(4>) with respect to (f>. The last term in Eq. V^, G], 
is the sum of all two-particle irreducible vacuum diagrams where all lines represent full propagators 
G. The expectation values of the one-point function, <p, and of the two-point function, Q{k) are 
determined from the stationarity conditions 



6V\<b,G\ 



5V[cj), G] 
8G(k) 



= 



,G=g 



= . 



With Eq. ©, the latter can be written in the form 

g-\k) = G- 1 (fc;0) + n(fc) 

where 



SG(k) 



(5a) 
(5b) 

(6a) 
(6b) 



!,=</,, G=g 



is the self-energy. Since n(fc) is in general a functional of Q, Eq. I|6a|l represents a Schwinger-Dyson 
equation for the full (dressed) propagator. 

The standard effective potential for the expectation value of the one-point function, V(4>), is 
obtained from the effective potential by taking the full propagator G to be a function of <fi, 
instead of an independent variable, 



where G(k; <p) is determined from 



V$) = V[4>,G<M 
8V[4>,G] 



, 



G=G 



which is equivalent to 



where 



SG{k) 

G- 1 {k-4>) = G- 1 {k-4>) + tl{k-< 
5V 2 [4>,G] 



n(fc;0) = 2 



5G{k) 



(7) 



(8) 



(9a) 



(9b) 



G=G 



This expression and Eq. @ can be used to obtain a compact form of the standard effective potential 



V$) = U{4>) + \ j In G" 1 (k; $) - \ j U(k; $) G(k; $) + V 2 [4>, G] . (10) 

Since V% contains infinitely many diagrams, an exact calculation is impossible. In practice, one has to 
restrict the computation of V% to a finite number of diagrams. The selected set of diagrams defines a 
particular many-body approximation. Cutting internal lines in these diagrams according to Eq. Ij6bfl , 
one obtains the diagrams contributing to the self-energy of the quasiparticles in this approximation 
scheme. Solving the Dyson-Schwinger equation H6a(l provides a self-consistent calculation of this 



G 




FIG. 1: Left-hand side: the double-bubble diagrams. Full lines are scalar particles, dashed lines are pseu- 
doscalar particles. Right-hand side: the tadpole contributions to the self-energies obtained by cutting a line 
in the double-bubble diagrams on the left-hand side. 



self-energy. In general, the Dyson-Schwinger equation is an integral equation for the self-energy as 
a function of energy and momentum. 

If one only takes the double-bubble diagrams on the left-hand side of Fig. ^ into account in the 
calculation of V%, one obtains the Hartree approximation. Cutting these diagrams yields the well- 
known tadpole diagrams for the self-energy, cf. the right-hand side of Fig. ^ The Dyson-Schwinger 
equation (|6a[) is a self-consistency equation for this self-energy due to the fact that the internal lines 
in the tadpole diagrams represent full propagators. The Hartree approximation is a particularly 
simple many-body approximation scheme, because the tadpole diagrams are independent of energy 
and momentum, and thus the Dyson-Schwinger equations are no longer integral equations, but 
become fix-point equations for the quasiparticle masses. 

B. The 0(4) model 

In this subsection, we apply the CJT formalism to the O(N) model. Our numerical results 
presented in Sec. IIVI are exclusively for the case N = 4. The CJT effective potential for the O(N) 
model is [Tcj 

V[*,S,P] = U(a) + ~ f [In S' 1 (fc) + S- 1 (k; a) S{k) - l] 

^ Jk 

+ J k [lnP- 1 (fc)+P- 1 (fc;a)P(fc)-l] + V 2 [a,S,P] , (11) 

where a is the expectation value of the scalar field (in the presence of external sources). Because 
the vacuum of QCD has even parity, the expectation values of the pseudoscalar fields can be set to 
zero. The quantities S and P are the full propagators for scalar and pseudoscalar particles, while 
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S 1 and P 1 are the corresponding inverse tree- level propagators, 



where the tree-level masses are 



ml(a) 



-k 2 +ml(a) , 
-k 2 + m 2 n {a) , 



o 12 A o 

+ ■ 



(12a) 
(12b) 



(13a) 
(13b) 



The constant /i 2 is the bare mass term in the Lagrangian of the 0(N) model, while A is the four-point 
coupling constant. For fi 2 < 0, the O(N) symmetry is spontaneously broken to 0(N — 1), leading 
to N — 1 Goldstone bosons. The tree-level potential is 



(14) 



where H is a term which breaks the O(N) symmetry explicitly to 0(N — 1). V^[<7,>S, P] is the 
sum of all two-particle irreducible diagrams. In the following, we restrict ourselves to the Hartree 
approximation, i.e., we take into account only the double-bubble diagrams shown on the left-hand 
side of Fig. ^ These diagrams have no explicit a dependence. Then, only tadpole diagrams (with 
resummed propagators) contribute to the self-energies, cf. the right-hand side of Fig. ^ In the 
Hartree approximation, 



V 2 [S,P] = 3 A 



S(k) 



+ (iV + l)(iV-l)A 



P(k) 



S(k) 



Hp) 



The stationarity conditions (|5a|l and l|5b|l read 

H = ^ a +^a 3 + ^aJ [3 S(q) + (N — 1) V(q)] 



V- X (k) = -k 2 



M 2 , 



(15) 



(16a) 

(16b) 
(16c) 



where M a and M n are the a and pion masses dressed by contributions from the diagrams of Fig.^ 



Ml 
Ml 



m 2 (cr) 
m 2 (cr) 



4A 
~N 
4A 
~N 



S(q) + (N-1) / P(q) 

Jq 

S(q) + (N+l) f V{q) 



Using Eqs. Ijl7a(l and l|17bj) . Eq. Ijl6a(l can be written in the compact form: 



H = a 



M 2 (T z 

TV 



(17a) 
(17b) 

(17c) 



Equations 117al) . (|17bl) . and H17c|l are the stationarity conditions of the O(N) model in the Hartree 
approximation. The explicit calculation of the tadpole integrals J S(q) and J V(q) will be discussed 
in Sec. dSl 



8 



C. The U{Nf) r x U(N f ) e linear sigma model for Nf = 2,3 and 4 flavors 



The application of the CJT formalism to the U(Nf) r x U(Nf)e linear sigma model for Nf = 3 
was discussed in Ref. [j| . Since we want to treat the cases Nj — 2 and Nj = 4 on the same footing, 
we derive the CJT effective potential in somewhat greater detail than in the last subsection. 

The Lagrangian of the U(Nf) r x U(Nf)t linear sigma model for Nf — 2, 3 or 4 flavors is given by 



(18) 



I '). 20. 21. 2 j 

£($) = Tr (d^d*® - to 2 $+$) - Ai [Tr ($ f $)] 2 - A 2 Tr ($t$ 

+ c [det ($) + det ($*)] + Tr [if ($ + . 

$ is a complex Afy x Nf matrix parametrizing the scalar and pseudoscalar mesons, 

— T a (f> a = T a K + m a ) , (19a) 

where a a are the scalar (J p = + ) fields and 7r a are the pseudoscalar (J p = - ) fields. The Nf x Nf 
matrix H breaks the symmetry explicitly and is chosen as 

H = T a h a , (19b) 

where h a are external fields. T a are the generators of U(Nf). The T a are normalized such that 
Tr(T a T b ) = S ab /2. They obey the U(N f ) algebra with 



[T a , Tb] 
{Ta,T b } 



i fabc T c j 



(20a) 
(20b) 



where / a f, c and <i a {, c are the standard antisymmetric and symmetric structure constants of SU(Nf), 
a,b,c= 1, 



,Nj - 1, and 



/a 




(20c) 



The terms in the first line of Eq. (fTHfl are invariant under U(Nf) r x U(Nf)t = U(Nf) v x U(Nf)A 
transformations. The determinant terms are invariant under SU(Nf) r x SU(Nf)e = SU(Nf)y x 
SU(Nf)A, but break the (7(1)^ symmetry explicitly. These terms arise from the U(1)a anomaly 
of the QCD vacuum. The last term in Eq. (|18fl breaks the axial and possibly the SU(Nf)y vector 
symmetries explicitly. 

In the following we discuss the three different cases Nf — 2, 3, and 4 in detail. The identification 
of the a a and 7r a fields with the physical scalar and pseudoscalar mesons is given in Appendix 1X1 
A non- vanishing vacuum expectation value for $, 



T a a a , 



(21) 



breaks the chiral symmetry spontaneously. (If the vacuum does not break parity, the fields n a cannot 
assume a non-vanishing vacuum expectation value.) Shifting the $ field by this expectation value, 
the Lagrangian can be rewritten as 

C = - [d^Uad^CJa + <9 M 7T a <9'V a - O a{m%) afeCTft - TT a (mp)ab^b) 



4 4 

g Fabcd &d - S(Nf,3) Qabc + g S(Nf, 4) Qabcd &d 

4 H abcd <T d + 3 S(N f , 3) G abc - 4 S(N f , 4) Q 

abed O d 



OaObCc 
Ka^bOc 



- 2 



H abcd -S(N f ,4)g 

abed 



OaOb^c^d 



1 

3 

U{a) , 



abed (a a cr b a c ad + Tra^b^c^d) 



(22) 
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where 5(n, m) = 8 nm is the Kronecker delta and the fields a a and 7r are the fluctuations around the 
expectation values a a . The latter are determined from the condition dU(a)/da a = 0. The tree-level 
potential is 



U(a) 



m 
~2 

1 r 
3 



3 S(N f ,2)Q ab + 8(N f ,3)Q abc a c 



&a&b 



Fabcd + 5{Nf : 4) Qabcd 



VaVbVcVd — h a (T a 



The coefficients T abcd , Gab, Gabc, Qabcd-, and Habcd are given by 

J\ abed 

r r 

SaoSbO — SalSbl — (5 a 2<5b2 — <5a3&>3 



-J {Sabred + 5 a dhc + S ac S bc [) + (d a bnd n cd + d a d n d n bc + d acn d nb d) 



Gab 
Gabc 
Gabcd 



3 9 

dabc — ^ (Saodobc + 5bod a Qc + Scodabo) + ^ ^OOO^aO^hO^cO 



(23) 

(24a) 
(24b) 
(24c) 



c 

16 



Sabred + Sad$bc + Sac^bd ~ {dabndncd + d a dnd n bc + d acn d nb d) + 16 S a0 S b oS c o5 d o 
4 (SaoSbO^cd + SaoficoSbd + S a oSd0^bc + SboScO^ad + 5bo5dO$ac + ^dO^cO^ab) 



7~La 



bed 



+ V8 (Saodbcd + Sbodcda + Scoddab + Sdodabc) 
^ab^cd ~t~ (,d abn d ncc [ -\- facnfnbd fbcnfyiad*) 



(24d) 
(24e) 



The tree-level masses, [?n|]ab and [mp] a h, are given by 

[ m s(°-)] ab = m 2 5 ab -6[S{N f ,2)Gab + 5{N f ,3)GabcVc}+4 [Fabcd + S(Nf,A)Gabcd] O c &d ,(25a) 
[m 2 P (a)] ab = m 2 5 ab + 6[5{N f ,2)Gab + 5(N f ,3)GabcVc}+4 [Habcd — 8(Nf, 4:)Gabcd) Vc&d -(25b) 

In general, these mass matrices are not diagonal. Consequently, the fields (er a , 7r a ) in the standard 
basis of U (Nf) generators are not mass eigenstates. Since the mass matrices are symmetric and real, 
diagonalization is achieved by an orthogonal transformation, 



(26a) 
(26b) 
(26c) 



\rh 2 1 - (s ' p) \m 2 1 (s ' p) 

\J n S,P\ l - U ai Y m S,P\ ab U bi 

The effective potential of the U(Nf) r x U(Nf)e linear sigma model in the CJT formalism reads 
V[a,S,P] = U(*) + lJ k ' {[\nS-\k)] aa +[]np-\k)] aa } 

+ 11 [Sab( k > °) Sba{k) + P- b \k- a) Pba(k) - 2S ab S ba ] + V 2 [a, S, P] . (27) 



(28a) 
(28b) 



Here, U(a) is the tree-level potential of Eq. Q23[) . and 

Sabi^a) = ~k 2 S ab +[m 2 (a)] ab , 
P*fc*) = -k 2 6 ab +[m 2 P (a)] ab , 



are the tree-level propagators for scalar and pseudoscalar particles, with the respective mass matrices 
(|25p. The fluctuations a a and 7r a around the expectation values a a no longer occur in the effective 
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potential (|27|) . Therefore, from now on we use the symbol a a for the expectation values for the 
scalar fields in the effective potential l|27|) . These expectation values, and the full propagators for 
scalar, S(fc), and pseudoscalar, V(k), particles are determined from the stationarity conditions 



5V[a, S, P] 



Scr a 
SV[a,S,P] 



a=a,S=S,P=V 



SSab 



o , 







6V[a,S,P] 



SPab 



g=cr,S=S,P=V 

With Eq. (|6b|l . the latter two equations can be written in the form 

= S- b \k-a) + ^ ab {k) , 

Kb( k ) = P^{k-a) + n ab {k) , 

where 

n SV 3 [a,S,P] 



g=a,S=S.P=V 



n ab (fc) 



5Sba(k) 

SV 2 [a,S, P] 



5P ba (k) 



g=a,S=S,P=V 



g=cr,S=S,P=V 



(29a) 
(29b) 



(30a) 
(30b) 



(31a) 
(31b) 



are the self-energies for the scalar and pseudoscalar particles. As in the case of the O(N) model, we 
include only the two- loop diagrams of Fig. ^ in V*. Then, 



V 2 [S,P] = [Tabcd + S(N f A)Ga 



bed 



Sab(k) / S cd ( P ) + / P ab (k) / P cd ( P ) 



+ 2[H a bcd-8{N f A)Gabcd\ / § ab {k) / P cd {p) . 



(32) 



Note that, in the Hartree approximation, V 2 is independent of a a . The stationarity conditions for 
the condensates are 

h a = m 2 a a - [6 S(N f , 2) Q ab + 3 S(Nf,3) Gabc <J C ] <?b + 7; [Fabcd + S(Nf,4) Gabcd] cr b a c a d 



(33) 



In the Hartree approximation, the self-energies l|6b() are independent of momentum, and the 
Schwinger-Dyson equations l|30a|) for the full propagators assume the simple form 

(34) 
(35) 



+ {-Z8{N f ^)g abc + A[T abc d + 5{NfA)Gabcd}od} / S cb (k) 

Jk 

+ {3 6(N f ,3)g abc + 4[H bcad -6(Nf,4)Gabcd\<Td} I ' V cb (k) . 



Sab(k) = -k 2 S ab +[Ml] ab 

Kb(k) = -k 2 5 ab +[M 2 P ] ab 
The scalar and pseudoscalar mass matrices are given by 



M 



Sl ab 



+ 4 [T abcd + S(N f , 4) G abcd ] / S cd (q) + 4 [H abc d - 5(Nf,4) Gabcd] / V cd {q) , (36a) 

J q J q 

[Ml] ab = [ml{a)] ah 

+ 4 [T ab cd + 5(Nf,4) Gabcd] I Vcd(q) + 4 [Habcd - S(Nf,4) Gabcd] I S cd (q) . (36b) 

J q J q 
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In the Hartree approximation, all particles are stable quasiparticles, i.e., the imaginary parts of the 
self-energies vanish. Therefore, the inverse propagators H34|) and (|35|1 are real-valued. They are also 
symmetric in the standard basis of U(Nf) generators and thus diagonalizable via an orthogonal 
transformation. This transformation is given by Eq. (|26c[l . with the obvious replacements 



KpU-> [ M lp] ab > [™%p] 



Mlp 



(37) 



The propagator matrices are diagonalized by the same orthogonal transformation as their inverse. 
The tadpole integrals in Eqs. (|33|l and l|3fci|> are therefore computed as 



s b M = o<£ ) I s t ( q )o 



V bc (q) = O 



-n(-P) 



Vi(q)o. 



(p) 



(38) 



After this rotation, only tadpole integrals over propagators with a single index have to be computed. 
This will be discussed in the next section. 



D. Explicit calculation of loop integrals and the effective potential 



In principle, the calculation of the tadpole integrals J^Si(q), J q S(q) and J^'Pi(q), and J V(q) 
requires renormalization. Renormalization of many-body approximation schemes is nontrivial, but 
does not change the results qualitatively 0, fiol ] . We therefore simply omit the vacuum contributions 
to the loop integrals, 



Si(q) 
V{q) 



d 3 q 1 
(2^F e q [(Mj)i] 

d 3 q 1 



exp ■ 



J.™ = / 



(2tt)3 £q [(M|) 
d 3 q 1 

d 3 q 1 



exp ■ 



£q[(Mg)j] 

T 



T 



- 1 



- 1 



(2tt)3 £q [M2] 



exp 



exp 



T 
T 



- 1 



1 



(39a) 

(39b) 
(39c) 
(39d) 



Here, e q [M 2 ] = -\/q 2 + M 2 is the relativistic energy of a quasiparticle with mass M and momentum 

q- 

Now we compute the standard effective potential V(a) = V[a, S(a), P(a)] from Eq. I|1U|) . Since 
V% has the general structure 



V 2 [S,P] =c s 



S(k) 



P(k) 



S(k) 



Uk 



Pip) 



cf. Eqs. I|15|l and (|27[1 . the self-energies l|31|) assume the form 

8Vi[<t,S, P] 



T,(k;a) = 2 
Il(k;a) = 2 



*Va [a, S, P] 



a=a,S=S,P=P 



\c s \ S(q;a) + 2c sp / P(q:a 



8P{k) 



a=a,S=S,P=P 



= Ac p / P(q;a) + 2. 



Sp 



S(q;a) 



With these expressions one derives the identity 



E(/c; a) S(k; a) 



1 



Il(k;a)P(k;a) = -2V 2 [S(a),P(a) 



(40) 

(41a) 
(41b) 

(42) 



12 



This considerably simplifies the expressions for the standard effective potential. For the O(N) model 
we obtain 

1 f N — 1 f - 

V(a) = U(a) + -J InS-^kia) + — ^— J \nP-\k-a) - V 2 [S(a), P(a)} , (43a) 

while for the U{Nf) r x U{Nf)t model we have 




(43b) 



The momentum integrals in Eqs. 143a|) and (|43bjl require renormalization, too. As above, we simply 
omit the vacuum contribution, which leads to the following integrals 0: 

Jk 

f \nPr\k-a) 

Jk 

J In S^ik; a) 
[ \nf-\k-d) 

Jk 

In these equations, the tilde denotes quantities which are diagonalized according to Eqs. H2(jcfl and 
(|38|l . respectively. A hat denotes propagators computed according to Eq. (|9a|) . Masses with a hat 
are the mass terms appearing in these propagators. 

III. PATTERNS OF SYMMETRY BREAKING AND VACUUM PROPERTIES 

In this section, we discuss the patterns of symmetry breaking in the vacuum and determine the 
coupling constants of the models from the vacuum properties of the mesons. 

A. Patterns of symmetry breaking 

For the 0{N) model (with N — 4) we study two different patterns of symmetry breaking (cf. 
Table |U): 

1. H = 0: For [x 2 < the O(N) symmetry is spontaneously broken to 0(N — 1), giving rise to a 
non-vanishing expectation value for the a field and N — 1 Goldstone bosons. 

2. H 7^ 0: The term —Ha in Eq. (|14|) corresponds to nonzero quark masses in the QCD La- 
grangian. It breaks the O(N) symmetry explicitly to 0(N — 1). The N — 1 Goldstone bosons 
become pseudo-Goldstone bosons. 

For the U(Nf) r x U(Nf)i models with Nj = 2 and 3 we study the following patterns of symmetry 
breaking (cf. Table QJ: 
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1. h a = 0,c = 0: For m 2 < the global SU(N f ) v x U(N f ) A symmetry is broken to SU(Nf) v , 
and <& develops a non- vanishing vacuum expectation value, ($) = Tq<7q. By the Vafa-Witten 
theorem [23J, only the axial symmetries can be spontaneously broken, while the vector symme- 
tries stay intact. In order to retain an SU(Nt)v symmetry, only the term proportional to To 
survives in the sum over a in Eq. (|21|) for the vacuum expectation value ($). Spontaneously 
breaking U(Nf)A leads to Nj Goldstone bosons which form a pseudoscalar, N'j dimensional 
multiplet. This case is referred to as the chiral limit without U(1)a anomaly. 

2. h a = 0, c ^ 0: The symmetry is SU(Nf)v X SU(Nf)A- A non-vanishing ($) spontaneously 
breaks the symmetry to SU(Nf)v, with the appearance of Nj — 1 Goldstone bosons which 
form a pseudoscalar, Nj — 1 dimensional multiplet. The Nj th pseudoscalar meson is no longer 
massless, because the U(1)a symmetry is already explicitly broken. This case is referred to as 
the chiral limit with U(1)a anomaly. 

3. h a =/= 0, c = 0: In QCD this corresponds to non-vanishing quark masses, but a vanishing U(1)a 
anomaly. Since ($) must carry the quantum numbers of the vacuum, only the fields a a 
corresponding to the diagonal generators of U(Nf) can be nonzero. The same holds for the 
external fields h a which generate a non-vanishing expectation value by explicitly breaking the 
U(Nj)a symmetry. For Nf = 2 these are h and h 3 , for Nt = 3 there is an additional field, h s . 
Because the masses of the up and down quarks are approximately equal, m u ~ ma, we restrict 
our study to ho ^ and hs = for all cases considered. Since the strange quark mass m s is 
larger than m u ~ irid, h% ^ 0. In this case the SU(Nf)v X U(Nf)A symmetry is explicitly 
broken to SU(2)v- The latter symmetry remains intact because /13 — 0. This case is referred 
to as the case of explicit chiral symmetry breaking without U(1)a anomaly. 

4. /i n ^, c ^ 0: A U(1)a subgroup of the U(Nf)A symmetry is explicitly broken by instantons. 
As explained above we restrict ourselves to /13 = 0. This case is referred to as the case of 
explicit chiral symmetry breaking with U(1)a anomaly. 

For Nf — 4, we only study the last case, i.e., explicit chiral symmetry breaking with U(1)a 
anomaly. To break the U (4)a symmetry explicitly, in addition to ho ^ and /is / we have 
to introduce a nonzero value for the field /115 corresponding to the fourth diagonal generator of 
C/(4). Since the charm quark mass is much larger than the light up and down quark masses and the 
strange quark mass, h 15 is also much larger than either h or h 8 , cf. Table Ivl Therefore, it does not 
make too much sense to study the rather unrealistic first two cases of vanishing quark masses. We 
therefore restrict our considerations to the physical case of explicit chiral symmetry breaking with 
U(1)a anomaly. 

B. Condensates and masses in the vacuum 

In this section we determine the parameters of the different models from the vacuum values of 
the condensates and the meson masses for the various symmetry breaking patterns discussed in Sec. 
IIII Al For the 0(4), the U(3) r x U(3)e, and the U(A) r x J7(4)^ linear sigma model, we simply follow 
Refs. 0, |22, respectively. Since it has not been done previously, fitting the parameters of the 
U(2) r x U(2)e model is discussed in more detail. 

For the 0(4) model, we have three parameters, H, A, and which are adjusted to reproduce the 
vacuum values for the pion decay constant, f w , the pion mass, m^, and the a mass, m a . For reasons 
explained below, for the latter we choose m a = 400 MeV, instead of 600 MeV as in Ref. 0. The 
values for the parameters in the chiral limit and with explicit symmetry breaking are listed in Table 

eh 

In the U(2) r x U{2)i model, for all symmetry breaking patterns studied here there is an (ap- 
proximate) SU(2)v symmetry due to the (approximate) equality of the up and down quark masses. 
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Masses and decay constants 


Parameter set 


Explicit chiral 

symmetry 

breaking 


f n = 92.4 MeV 
m CT = 400 MeV 
m n = 139.5 MeV 


H = (121.60 MeV) 3 
A = 8.230 

H 2 = -(225.41 MeV) 2 


Chiral limit 


,U = 90 MeV 
m CT = 400 MeV 
m n = MeV 


H = 
A = 9.877 

y 2 = -(282.84 MeV) 2 



TABLE II: The masses and decay constants at vanishing temperature and the corresponding parameter 
sets for the 0(4) linear sigma model for the two symmetry breaking patterns studied here. 



Consequently, for all cases the vacuum expectation value is ($) = ToCTo- At zero temperature, the 
equation for the condensate (Tq reads, cf. Eq. . 



ho = (To 



m 2 -c+(\ 1 + ^) a 2 



(45) 



The PCAC relations determine the value of the condensate from the pseudoscalar meson decay 
constants, 

fa = daaO C ■ (46) 

Since d aa o = lj ah meson decay constants are identical, and we obtain ao = fw The scalar mass 
matrix is diagonal, with the elements 



ml = 



m l 



[ m U^o)] 00 = m 2 ^c + S^M + ^a 2 , (47a) 
[m|(<7 )] n = [m|(£r )] 2a = [m|(<r )] 33 = m 2 + c + ( \ t + ^ ) a 2 . (47b) 



The pseudoscalar mass matrix is also diagonal, with the elements 



m 



2 ~- Kp(^o)] 00 = m 2 + c + (Ax + a 2 , (48a) 



>M U = [m|(ao)] 22 = [m 2 3 (a )] 33 = m 2 - C + (^Ai + ^^a 2 . (48b) 



Without the U(1)a anomaly, c = 0, the pions and the rj meson become degenerate in mass. In the 
chiral limit, one then has four (instead of three) Goldstone bosons. With the U(1)a anomaly, c is 
positive, cf. Table Mil and the r\ meson becomes heavier than the pion. At zero temperature, the 
(squared) mass difference between the rj and the pion is determined by the parameter c characterizing 
the strength of the U(1)a anomaly, m 2 — m 2 — 2c. Simultaneously, also the mass difference between 
the ao and the a meson is determined by this parameter, m 2 — m 2 = 2c — 2Ai<7 2 . (As Ai < 0, cf. 
Table nTTl the second term always increases the mass difference.) 

The limit c — * oo corresponds to maximum explicit U{1)a symmetry breaking. In this limit, for 
realistic values of the a meson and the pion mass (i.e., m 2 — c = const.), the n and ao mesons 
become infinitely heavy and are thus removed from the spectrum of physical excitations. In this 
limit, the U(2) r x U(2)i is identical to the 0(4) model, where the ao and n meson are absent from 
the beginning. 
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Masses and decay constants 


Parameter set 


Explicit chiral 
symmetry breaking 
with U(1)a 
anomaly 


f„ = 92.4 MeV 
ra a = 400 MeV 
m ao = 984.7 MeV 
= 139.5 MeV 
m„ = 547 MeV 


h = (121.60 MeV) 3 

Ai = -31.03 

A 2 = 78.52 

m 2 = (298.44 MeV) 2 

c = (374.00 MeV) 2 


explicit chiral 
symmetry breaking 
without U(1)a 
anomaly 


/„ = 92.4 MeV 

m CT = 400 MeV 

m ao = 984.7 MeV 

m n = m v = 139.5 MeV 


h = (121.60 MeV) 3 

Ai = -47 '.41 

A 2 = 111.29 

m 2 = -(225.41 MeV) 2 

c = 


Chiral limit 
with U(1)a 
anomaly 


U = 90 MeV 
m a = 400 MeV 
m ao = 984.7 MeV 
m w — 

m, = 547 MeV 


h =0 

Ai = -31.51 

A 2 = 82.77 

m 2 = (263.83 MeV) 2 

c= (386.79 MeV) 2 


Chiral limit 
without (7(1) a 
anomaly 


U = 90 MeV 
m a = 400 MeV 
m ao = 984.7 MeV 
m w — rriri = 


ho = c = 

Ai = -49.98 

A 2 = 119.71 

m 2 = -(282.84MeV) 2 



TABLE III: The masses and decay constants at vanishing temperature and the corresponding parameter 
sets for the U(2) r x U(2)t model for the four symmetry breaking patterns studied here. 



With Eqs. 145|) . I|47|l . and (|48|l . we can determine the parameters of the model from the pion decay 
constant and the meson masses in the vacuum, 



CO — fir i M 



11 11 

„ 2 — m% m — 

m = m„ 4 ! — — 



2/| ' fl ' 

2 2 

c = m "~ m " , = Uml . (49) 



With the U(1)a anomaly, there are five parameters, ho, Ai, A2, m 2 , and c, which can be unambigu- 
ously determined from the five quantities mo-, m ao , m v , and Without the U(1)a anomaly, 
c = 0, and m v = m w . In this case, there are only four parameters and four quantities from which 
the values of the parameters can be fixed. The values for the parameters are listed in Table ITTT1 

For the U (3) r x U (3)t model, we follow Ref. in fitting theparameters of the model to vacuum 
quantities. Our parameters differ from the ones given in Ref. [9|, since we use m a = 400 MeV, and 
not 600 MeV. In the chiral limit, the number of parameters equals the number of vacuum quantities, 
and one can again obtain a unique mapping between these sets of quantities. With explicit chiral 
symmetry breaking, however, there are fewer parameters than vacuum quantities. Consequently, 
some meson masses are predicted rather than used as fit parameters. The values for the parameters 
and the meson masses are given in Table HVl The vacuum quantities predicted by the fit are given 
in bold-faced letters. 

For the U(4) r x U(4)e model we adjust the parameters to obtain reasonable agreement between 
the vacuum masses [24| and the masses computed at tree-level, and not the masses computed to 
one-loop order as in Ref. j^. As for the f/(3) r x U(3)e case, the number of parameters is smaller 
than the number of meson masses, such that some meson masses cannot be fitted independently, 
but are predicted within this approach. We found that small values for the a meson, m a ~ 400 
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Masses and decay constants 


Parameter set 


Explicit chiral 
symmetry beaking 
with U(V) A 
anomaly 


U = 92.4 MeV 
f K = 113 MeV 
m CT = 400 MeV 
m ao = 1024.6 MeV 
m K = 1116.2 MeV 
m /o = 1188.7 MeV 

= 139.5 MeV 
m K = 493 MeV 
m, = 536.5 MeV 
m,. = 963.9 MeV 


ho = (285.04 MeV) 3 

h 8 = -(309.46 MeV) 3 

Ai = -5.38 

A 2 = 45.08 

m 2 = (493.69 MeV) 2 

c = 4831.25 MeV 


Explicit chiral 
symmetry breaking 
without U(1)a 
anomaly 


f v = 92 A MeV 
f K = 113 MeV 
m a = 400 MeV 
m ao = 844.4 MeV 
m K = 1116.2 MeV 
m /o = 1248.3 MeV 
m w — m v i — 139.5 MeV 
m K = 493 MeV 
m, = 630.6 MeV 


ho = (285.04 MeV) 3 

hg = -(309.46 MeV) 3 

Ai = -24.13 

A 2 = 81.24 

m 2 = (306.50 MeV) 2 

c = 


Chiral limit 
with U(l) A 
anomaly 


U = /a- = 90 MeV 
m„ = 400 MeV 

m ao = m K = m/ = 1225.8 MeV 
m-n- — m,K ~ m v — 
m v > = 958 MeV 


ho = h 8 = 

Ai = -17.48 

A 2 = 109.97 

m 2 = (270.11 MeV) 2 

c = 6798.25 MeV 


Chiral limit 
without U(1)a 
anomaly 


U = /A- = 90MeV 
m CT = 400 MeV 

m ao = m K = m/ = 1225.8 MeV 
m,r = m,K ~ m-q = m„/ = 


ho = hs — c = 

Ai = -55.25 

A 2 = 185.50 

m 2 = -(282.84 MeV) 2 



TABLE IV: The masses and decay constants at vanishing temperature and the corresponding parameter sets 
for the U(3) r x U(3)e model for the four symmetry breaking patterns studied here. The masses in bold-faced 
letters are predicted, the other masses and decay constants are used to calculate the parameter set. 

MeV, are favored, otherwise the mass spectrum of the charmed mesons deviates too much from the 
one in nature. This is the reason why we choose a a meson mass m a — 400 MeV also in the other 
cases discussed above. The values for the parameters and vacuum quantities are listed in Table IVI 

IV. RESULTS 

In this section, we discuss the numerical results at nonzero temperature for the cases listed in 
Table H| 
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FIG. 2: The meson masses as a function of temperature for the different models studied here, for the case 
with U(1)a anomaly and explicit chiral symmetry breaking. 
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Masses and decay constants 


Parameter set 


Explicit chiral 


f v = 92 .4 MeV 


ho = 


(917.24 MeV) 3 


symmetry breaking 


f K = 113 MeV 


h s = 


-(309.46 MeV) 3 


with U(1)a 


m CT = 400.6 MeV 


his = 


= -(1088.67 MeV) 3 


anomaly 


m ao = 1052.6 MeV 


Ai = 


-0.12 




m K = 1116.2 MeV 


A 2 = 


4.85 




m /a = 1178.6 MeV 


m 2 = 


= (345.78 MeV) 2 




m Do = 2370.0 MeV 


c = - 


-1.50 




m Ds0 = 2480.1 MeV 








m Xc0 = 3565.6 MeV 








m„ = 139.5 MeV 








m K = 493 MeV 








m, = 542.5 MeV 








m n , = 1028 MeV 








m D = 1944.1 MeV 








m Ds = 1899.1 MeV 








m^ c = 2129.8 MeV 







TABLE V: The masses and decay constants at vanishing temperature and the corresponding parameter sets 
for the U(4) r x U{A)i model for the case of explicit chiral symmetry breaking with U(1)a anomaly. The 
masses written in bold-faced letters are predicted, the other masses and decay constants are used to calculate 
the parameter set. 

A. Explicit chiral symmetry breaking with U(1)a anomaly 

In Fig. [2 we show the masses of the mesons as a function of temperature for explicit chiral 
symmetry breaking, including the U(1)a anomaly. This is the case where chiral symmetry breaking 
results in the smallest residual symmetry group, SU(Nf)v x U(Nj)a — > SU(2)v- 

In Fig. |H^a) the masses of the a meson and the pions are shown for all models. The a meson 
and the pion become degenerate in mass in the chirally restored phase. Comparing the a meson 
and pion masses in the 0(4) model with those in the U(2) r x U{2)i model, the difference is almost 
negligible up to temperatures of 150 MeV. In the chirally restored phase, the masses behave linearly 
with temperature, but grow faster in the U(2) r x U(2)t model than in the 0(4) model. The reason 
is that there are twice as many fields in the former model than in the latter, which results in twice 
as many tadpole-like contributions in the equations for the in-medium masses. These come with a 
positive sign and thus increase the masses. 

Comparing the results of the U(3) r x U(3)t model with those of the U(2) r x U(2)g model, one 
observes differences already at a temperature of about 100 MeV. Furthermore, the masses become 
even larger in the chirally restored phase. The reason for this behavior are the strange degrees of 
freedom in the Nf — 3 case which lead to additional tadpole-like terms in the self-energies. As 
above, they lead to an increase of the in-medium masses. 

Finally, one observes that virtually nothing changes in the temperature range of interest when 
including the charm degrees of freedom in the framework of the f/(4) r x f/(4)^ model. The reason 
is that the charm quark is large compared to the temperature, m c 3> T, and the contributions 
from charmed particles to the equations for the in-medium masses is suppressed. For two reasons, 
this is a non-trivial result. First, the equations for the in-medium masses are structurally different 
for the U(4) r x U(4)t model as compared to the U(3) r x U(3)t model, cf. Eqs. and 
Second, although the tadpole terms (|39|l arc strongly suppressed for particles with masses much 
larger than the temperature, Eqs. l\.V.'A) and 136|) form a nonlinear system of coupled equations, i.e., 
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small perturbations could lead to large quantitative changes in the solution. 

In Fig. Hfb) the masses of the ao and the 77 mesons are shown as functions of temperature. 
Qualitatively, the behavior of these masses is the same in all models. The eto meson mass is constant 
up to temperatures of 150 - 200 MeV. It then decreases, before increasing again above temperatures 
of 200 - 250 MeV. The r\ meson mass is constant up to T ~ 50 MeV, and then monotonously 
increases with temperature. At large temperatures, ao and 77 become degenerate in mass, indicating 
restoration of chiral symmetry. In the U{2) r x U(2)g model, this happens somewhat earlier, at about 
250 MeV, than in the other two cases. 

In the U(2) r x U(2)t model, the ao and 77 meson masses are used to determine the parameters 
of the model. Thus, at zero temperature, the masses coincide with their correct vacuum values. In 
the U(3) r x U(3)e model, the predicted 77 mass deviates only by 2% from its vacuum value. The 
ao mass is also rather close to the correct value; the predicted mass is about 4% too large. In the 
U(4) r x U(4)e, the 77 mass is reproduced with excellent accuracy (the deviation is less than 1%), 
while the ao mass is within 7% of its vacuum value. 

In Fig. |2{c) the K meson (now referred to as Xq(1430) 24]) and kaon masses are shown as a 
function of temperature. The results for the £7(3) r x U(3)i and J7(4) r x U(4)t model are almost 
identical. The k meson and kaon become degenerate in mass at temperatures of the order of 400 
MeV. In both models, the vacuum kaon mass is used as input, while the k mass is predicted. The 
deviation to the vacuum value in nature is about 21%. 

Figure [2Jd) shows the masses of the /o and the 77' meson as a function of temperature. These 
mesons also become degenerate in mass at temperatures of the order of 400 MeV. In the U(3) r x U(3)e 
model, the predicted rj 1 mass is rather close to its value in nature; the deviation is 0.6%. In the 
U(4) r x U(4)i model, the rf mass deviates from its correct vacuum value by about 7%. The /o 
mass is predicted in both models. If we identify this state with the /o(1370), these predicted masses 
deviate by 14% from their correct values. 

The masses of the D s q, Do, D Sl and D mesons are shown in Fig. Pfe) . The masses of the pseu- 
doscalar mesons are known but the scalar mesons and their masses have not yet been experimentally 
identified. It is somewhat peculiar that the charmed, strange D s meson is lighter than the charmed, 
non-strange D meson. This is an artifact of the particular set of coupling constants chosen here. For 
a different choice, this unphysical ordering of the masses can be reversed. Then, however, the masses 
of the other mesons deviate by an unacceptably large extent from their physical values. Temperature 
has virtually no effect on the heavy charmed mesons: their mass changes at most by 10%, even in 
the chirally restored phase. Due to the non-linear nature of the coupled system of Eqs. I|33|l and 
(I3(j|) . this is a non-trivial result, although not completely unexpected: we expect significant changes 
of the meson masses only when the temperature becomes of the order of the mass. For the charmed 
mesons with masses of the order of 2 GeV, this is never the case in the temperature range of interest. 

Finally, in Fig.[2{f) we show the masses of the Xco and ?7 C meson. Their large masses do not change 
at all for temperatures below 450 MeV. The vacuum values for both meson masses are predicted. 
While the mass for the Xco is within 5% of its correct value, the deviation for the r) c is somewhat 
larger (~ 30%). 

To summarize, the masses of the scalar mesons remain approximately constant up to temperatures 
around 150 MeV, and then slightly decrease before they become degenerate with the masses of the 
pseudoscalar mesons. On the other hand, with the exception of the 77' mass, the pseudoscalar masses 
in general increase monotonously with temperature. 

In Fig. EUa) the up-down quark condensate and (b) the strange quark condensate are shown as 
functions of temperature. In the 0(4) model and in the U(2) r x U(2)t model the up-down quark 
condensate y U p-down can be directly identified with the vacuum expectation value of the a field, 
^up— down = cr . On the other hand, in the U(3) r x U(3)t model 
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FIG. 3: (a) The up-down quark condensate and (b) the strange quark condensate as functions of temperature 
for the different models in the case with U(1)a anomaly and explicitly broken chiral symmetry. 



and in the U (4) r x U(4)e model the condensates are given by 
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In these formulas Estrange and y c harm are the strange and charm quark condensate, respectively. 

All models predict a qualitatively similar behavior for the temperature dependence of the up-down 
quark condensates. The strange quark condensate decreases more slowly with temperature than the 
up-down quark condensate. This is what one intuitively expects, as it appears more difficult to 
"melt" a condensate of heavier quark species than of light quark species. 

Figure 0] shows the charm quark condensate. Note that this condensate is much larger than 
the other two condensates. A peculiar feature is that it first increases at a temperature of about 
200 MeV, assumes a maximum at about 400 MeV, and then decreases. The maximum value is 
approximately 10% larger than the vacuum value. We are not aware of a simple explanation for this 
behavior. 



B. Explicit chiral symmetry breaking without U(1)a anomaly 

In Fig. [SI we show the masses for the scalar and pseudoscalar mesons for the case of explicit chiral 
symmetry breaking in the absence of the U(1)a anomaly. We discuss the results in comparison to 
the previous case. As in the previous case, the scalar meson masses stay constant up to temperatures 
close to the transition region, then decrease and finally start to increase again when they become 
degenerate with the pseudoscalar masses. In general, the pseudoscalar masses increase monotonously 
with temperature. The difference between the results obtained in the U(2) r x U(2)i and U(3) r x U (3)i 
model is rather small. An exception is the ao mass which is an input parameter in the former model, 
but is predicted in the latter. 
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FIG. 4: The charm quark condensate as a function of temperature. 



A marked difference to the case with U{\)a anomaly is that the chiral symmetry restoration 
transition is much more rapid, and it occurs at a slightly smaller temperature, T ~ 180 MeV. 
Moreover, above the transition the scalar and pseudoscalar masses become degenerate much more 
rapidly. The reason is the absence of explicit U(1)a symmetry breaking. At small temperatures and 
above the transition, the masses of the pion and the r\ meson are the same, because of the absence 
of the E/(1)a anomaly. In the temperature range from about 50 MeV to 225 MeV, however, they 
are different. We believe this to be an artifact originating from the violation of Goldstone's theorem 
in the Hartree approximation, which becomes even more obvious when considering the chiral limit. 

The melting of the condensates is shown in Fig. The smaller transition temperature of about 
180 MeV is also apparent in the temperature dependence of the up-down quark condensate. Again, 
the strange quark condensate melts less rapidly than the up-down quark condensate. 



The masses as a function of temperature for the 0(4) model are shown in Fig. Efa,l), for the 
U(2) r x U(2)t model in Fig. Eft),!.), and for the U(3) r x 17(3)* model in Fig. dc,!.). In the (9(4) 
and U(2) r x U(2)i models, there are three Goldstone bosons, the pions, while in the U(3) r x [7(3)* 
model there are eight Goldstone bosons, the pions, the kaons, and the r\ meson. The rj meson is not 
a Goldstone boson due to the explicit breaking of the U(1)a symmetry by the anomaly. The scalar 
octet, comprising the three ao mesons, the four k mesons, and the fo meson, is degenerate in mass, 
while the mass of the singlet a differs from the mass of the octet. As the temperature increases, the 
scalar masses decrease while the pseudoscalar masses increase. The mass of the Goldstone bosons 
increases, because the Hartree approximation does not respect Goldstone's theorem at nonzero 
temperature |^.ITo|. 

Due to the restoration of chiral symmetry above the transition temperature T c , the masses of the 
chiral partners become degenerate for temperatures T > T c . For the (9(4) model, the chiral partners 
are the a and the pion, for the U(2) r x U(2)e model they are the a and the pion, as well as the 
ao and the rj. Due to the explicit breaking of the U(1)a symmetry, the tr/pion and ao/rj do not 
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FIG. 5: The meson masses as a function of temperature for the different models in the case without U(1)a 
anomaly and explicit chiral symmetry breaking. 



become degenerate. The reason for this behavior is the term ~ Q a b in Eqs. Ij25(l . As discussed in Sec. 
IIII Bl at zero temperature this term leads to a difference in the masses (squared) of r\ meson and 
pion, and of the ao and a meson, respectively, which is proportional to 2c. In the case with U(\)a 
anomaly, 2c ^ even for temperatures above T c . Consequently, these mass differences persist also 
in the chirally restored phase. 

For Nf = 3, the situation is different. The term ~ Q a b for Nf = 2 is replaced by a term ~ Gabc&c- 
In the chirally symmetric phase, a c — 0, and consequently this term vanishes. Therefore, all meson 
masses become degenerate, even in the case with U(1)a anomaly. 

All models exhibit a first order phase transition between the low-temperature phase where chiral 
symmetry is broken and the high-temperature phase where chiral symmetry is restored. For the 0(4) 
model and the U(2) r x U(2)t model with explicit breaking of the U(1)a symmetry, the transition 
should be second order, cf. Sec.[I] It is a known shortcoming of the Hartree approximation to predict 
a first order transition also in these cases. For the U(2) r x U(2)g model without explicit breaking 
of the U(1)a symmetry, and for all U(Nf) r x U(Nf)g models with Nf > 2 the transition is of first 
order, which is correctly reproduced by the Hartree approximation. 

We have determined the numerical value of T c by computing the effective potentials. The latter are 
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FIG. 6: (a) The up-down quark condensate and (b) the strange quark condensate as a function of temperature 
for the different models in the case without U(1)a anomaly and explicitly broken chiral symmetry. 



shown for the three different models in the second column of Fig. [3 as a function of the condensate 
(Tq. (In the chiral limit, a is the only non-trivial condensate.) For the extraction of the transition 
temperature, the absolute normalization of the effective potential is irrelevant. All that matters is 
to identify the temperature where the minimum at the origin and the one at a nonzero value of <7rj 
become degenerate. For this purpose we have plotted the effective potential for two temperatures, 
one slightly below and one slightly above T c . From this we deduce that, for the 0(4) model, Fig. 
E[a,2), T c is between 159 and 160 MeV. For the U{2) r x U{2) e model, Fig. |7[b,2), we obtain a 
critical temperature between 154 and 155 MeV. This temperature is rather close to the one in the 
0(4) model. Finally, for the U(3) r x U(3)e model, the critical temperature is between 165 and 166 
MeV, which is slightly larger than in the previous cases. 

These values are surprisingly close to those obtained from lattice QCD calculations Q. In the 
chiral limit, these calculations predict T c ~ 175 MeV for N f = 2 and T c ~ 155 MeV for N f = 3. The 
critical temperature obtained from the chiral models in the Hartree approximation deviates from 
these values only by about 20 MeV (or 12%) for N f = 2 and only 10 MeV (or 6%) for N f = 3. 
However, for the chiral models the transition temperature is larger in the three-flavor case than in 
the two-flavor case, while one finds the opposite behavior in the lattice QCD calculations. 

For the sake of completeness we also show the condensates in the third column of Fig. [7\ Note 
that, in contrast to the cases where chiral symmetry is explicitly broken, the strange condensate 
is smaller than the up-down condensate. The reason is that, in the chiral limit, eg = 0, such that 
^Pup— down — V2 Estrange i cf. Eq. HJ. 



D. Chiral limit without U(1)a anomaly 

The masses as a function of temperature are shown in Fig. |H^a,l) for the U(2) r x U(2)e model, 
and in Fig.[Sfb,l) for the U(3) r x U(3)e model. In the first case, there are four Goldstone bosons, 
the pions and the rj meson, while in the latter case there are nine Goldstone bosons, the pions, the 
kaons, and the rj and 77' mesons. As the temperature increases, the scalar masses decrease while 
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FIG. 7: The chiral limit with U(1)a anomaly. The temperature dependence of the meson masses are shown 
in panels ( , 1) and that of the up-down and strange quark condensates in panels ( , 3). The results for the 
0(4) model are shown in panels (a, ), for the U(2) r x U(2)i model in panels (b, ), and for the U(3) r X U(3)i 
model in panels (c, ) . The effective potential (in arbitrary units) as a function of the condensate cro is shown 
in panels ( , 2). 

the pseudoscalar masses increase, until they become degenerate in a first order phase transition. 
The mass of the Goldstone bosons increases with temperature, because the Hartree approximation 
does not respect Goldstone's theorem at nonzero temperature [jj^j. Moreover, the masses of the 
Goldstone bosons are not equal: for the U(2) r x U(2)g model, the ij meson becomes heavier than 
the pion, while for the U(3) r x U(3)e model the rj meson becomes heavier than the other Goldstone 
bosons (pions, kaons, rj meson). Note that the role of the ij meson in the two-flavor case is assumed 
by the 77' meson in the three-flavor case. The reason is that the physical meson corresponding 
to the singlet representation in the U(2) r x U(2)t model is the ij meson, while it is the rj' in the 
U(3) r x 17(3)/ model. 

The numerical values for the critical temperature T c have been determined by computing the 
effective potentials. The latter are shown in the second column of Fig. Has a a function of the 



25 




50 100 150 200 250 300 350 400 20 40 60 80 100 120 50 100 150 200 250 300 



T(MeV) ct (MeV) T (MeV) 



FIG. 8: The chiral limit without U(1)a anomaly. The temperature dependence of the meson masses are 
shown in panels ( , 1) and that of the up-down and strange quark condensates in panels ( , 3). The results 
for the U(2) r x U(2)e model in panels (a, ), and for the (7(3) r x U(3)e model in panels (6, ). The effective 
potential (in arbitrary units) as a function of the condensate <to is shown in panels ( , 2). 

condensate ctq, for temperatures slightly below and above T c . In the U(2) r x U(2)e model, Fig. 
IHa,2) we obtain a critical temperature of ~ 150 MeV. For the U(3) r x U(3)i model the temperature 
is slightly smaller, between 147 and 148 MeV. The ordering of critical temperatures for the two- and 
three-flavor cases is now reversed as compared to the case with U(1)a anomaly, and consequently in 
agreement with the ordering found in lattice QCD calculations. One might be tempted to view this 
as an indication for a rapidly decreasing U(1)a anomaly near the critical temperature, in agreement 
with the results of Ref. [2ol f26j . Finally, the condensates are shown in the third column of Fig. [5] 
The results are qualitatively similar to the case with U(1)a anomaly. 

V. CONCLUSIONS 

In this work, we have used several different chiral models, the 0(4), the U(2) r x U(2)i, the 
U(3) r x U(3)e, and the U(4) r x U(4)t linear sigma model, to compute the temperature dependence 
of meson masses and quark condensates across the chiral phase transition. The meson masses and 
condensates were self-consistently calculated in the Hartree approximation, which we derived via the 
CJT formalism. Moreover, we have studied several distinct patterns of symmetry breaking within 
the different models. For a list of cases studied here see Table [|] 

We first considered the physically relevant case of explicit symmetry breaking in the presence of 
the U{\)a anomaly and compared the results of the different chiral models in order to clarify how 
they change with the number of quark flavors Nf. Comparing the 0(4) model with the U(2) r x U(2)i 



2G 



model, one first notices that the degrees of freedom have doubled: in addition to the a meson and 
the pions which are already present in the 0(4) model, one now has in addition the r\ meson and the 
ao mesons. This has the consequence that the meson masses grow more rapidly with temperature in 
the phase where chiral symmetry is restored. The reason for this are the tadpole contributions from 
the additional degrees of freedom to the meson self-energies, which lead to an increase in the meson 
masses. This result also applies when adding the strange degree of freedom in the framework of the 
U(3) r x U(3)i model. In fact, this picture holds in general, as long as the masses of the additional 
degrees of freedom are of the same order of magnitude as the chiral phase transition temperature. 
On the other hand, adding the heavy charm quark degree of freedom in the framework of the 
U(4) r x U(A)i model does not significantly influence the results for the masses of the non-charmed 
mesons and the non-charmed condensates. The reason is that the additional tadpole contributions 
from the heavy charmed mesons are exponentially suppressed with the meson mass, ~ exp(— M/T). 
Vice versa, also the masses of the charmed mesons do not change appreciably from their vacuum 
values over the range of temperatures of interest for chiral symmetry restoration, simply because the 
tadpole contributions from the non-charmed meson are small compared to the large vacuum mass of 
the charmed mesons. This result is intuitively clear from the physical point of view, but is still non- 
trivial: first, the equations for the in-medium masses are structurally different for the U(4) r x U(4)e 
model as compared to the U(3) r x U(3)e model. Second, the set of coupled equations for the masses 
and condensates is a nonlinear system of equations, which means that small perturbations could 
lead to large quantitative changes in the solution. 

We then studied the case of explicit chiral symmetry breaking without U{1)a anomaly. The main 
difference to the previous case was that the region of the chiral transition is narrower and located 
at a somewhat smaller temperature. 

Finally, we considered the meson masses and quark condensates in the chiral limit. The Hartree 
approximation correctly predicts the chiral transition to be of first order in the U{2) r x U{2)i model 
without U{1)a anomaly and in the J7(3) r x U(2>)i model. For the 0(4) model and the U(2) r x U(2)i 
model with U(1)a anomaly the Hartree approximation incorrectly produces a first order instead 
of a second order phase transition. The transition temperatures are surprisingly close to the ones 
obtained in lattice QCD calculations. However, in the case with U(1)a anomaly the transition 
temperature increases with the number of flavors, while in lattice QCD it decreases. This picture 
changes in the case with U(1)a anomaly, where the transition temperature shows the same behavior 
with the number of quark flavors as in lattice QCD. This may indicate that the U{1)a symmetry is, 
at least partially, restored at and above the chiral phase transition temperature. 
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For Nf = 2, the identification of the physical scalar and pseudoscalar meson fields with the matrix 
fields defined in Eq. I|19afl is 
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APPENDIX A: SCALAR AND PSEUDOSCALAR MESON FIELDS 




(Alb) 



(Ala) 



27 



Here, ^ = (tti =f i^)/^/^ and n° = 1T3 are the charged and neutral pions, respectively. Note the 
change of sign in the definition of the charged pion fields in terms of -k\j, in comparison to Ref. 

The definition given here is the correct one, as one can readily confirm by writing the meson 
fields in terms of their quark content, 



7T a ~ q T a 75 q 



(A2) 



This applies also to the other charged meson fields defined in the following. The field rj = ttq 
can be identified with the 77 meson. The parity partner of the pion is the ao(980) meson, i.e., 



+ _ 



(a 1 =F *o'2)/v / 2 and aj] = 03. The field ctq corresponds to the a meson [now also referred to as 



fo (400 - 1200)]. 

For Nf = 3 we obtain the following matrix: 
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The fields K ± = (tt 4 T « tt 5 )/\/2, if = (tt 6 - i tt 7 )/\/2, and if = ( 



)/v2 are the kaons. In 



7T6 + « 7T7 

general, because the strange quark is much heavier than the up or down quarks, the ttq and the ir$ 
are admixtures of the r\ and the rf meson. We identify the parity partner of the kaon with the k 
meson [now referred to as -£Tq(1430) in [24[]. Finally, in general the er and the <7g are admixtures of 
the a and Jo (980) mesons. 

For Nf = 4 the following identification of physical fields with matrix elements holds: 
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with 
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Here, £>° = {ir 9 + iTr w )/y/2, 5° = (7r 9 -Mn )/\/2, = (tth ±mti 2 )/n/2, and L>± = (^±^14)/^ 
are the charged and neutral pseudoscalar mesons with charm quantum numbers. The ttq, irg, and 
7Ti5 fields are admixtures of the 77, 77', and rj c mesons. The charged and neutral scalar mesons with 
charm quantum numbers arc D® = (er 9 + ia w )/V^, Dq = (erg — ia w )/V2, Dq = (an ± ia\2)/V^, 
and Df = (013 ± iai4)/v / 2- These scalar mesons have not been identified experimentally, yet. The 
(7 , <t$, and CT15 fields are admixtures of the a, fo, and Xco mesons. 
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